Signature of wave localisation in the time dependence of a reflected pulse 
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The average power spectrum of a pulse reflected by a disordered medium embedded in an iV-mode 
waveguide decays in time with a power law t~ p . We show that the exponent p increases from | to 
2 after N 2 scattering times, due to the onset of localisation. We compare two methods to arrive at 
this result. The first method involves the analytic continuation to imaginary absorption rate of a 
static scattering problem. The second method involves the solution of a Fokker-Planck equation for 
the frequency dependent reflection matrix, by means of a mapping onto a problem in non-Hermitian 
quantum mechanics. 
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The time-dependent amplitude of a wave pulse re- 
flected by an inhomogeneous medium consists of rapid 
oscillations with a slowly decaying envelope. The power 
spectrum a(u),t) describes the decay with time t of the 
envelope of the oscillations with frequency to. It is a basic 
dynamical observable in optics, acoustics, and seismology 
In the seismological context, the attention has fo- 
cused on randomly layered media, which are a model for 
the subsurface of the Earth. The fundamental result of 
White, Sheng, Zhang, and Papanicolaou H for this prob- 
lem is that a(u>, t) decays as t~ 2 for times long compared 
to the scattering time r s at frequency u>. The dynamics 
on this time scale is governed by localisation, since the 
product of t s and the wave velocity c equals the localisa- 
tion length in one dimension. Although this result for the 
power spectrum is more than a decade old, it has thus far 
resisted an extension beyond one-dimensional scattering. 

Work towards such an extension by Papanicolaou and 
coworkers HH has concentrated on locally layered me- 
dia, in which the scattering is one-dimensional on short 
length scales and three-dimensional on long length scales. 
This is most relevant for seismological applications. Re- 
cent dynamical microwave experiments by Genack et al. 
ra have motivated us to look at this problem in a wave- 
guide geometry, in which the scattering is fully three- 
dimensional — but restricted to a finite number N of 
propagating waveguide modes. (The single-mode case 
./V = 1 is statistically equivalent to the one-dimensional 
model of Ref. gj.) We find that the long-time decay 
of the average power spectrum is a power law as in the 
one-dimensional case, but with two exponents: a decay 
cx i~ 3 / 2 crosses over to a t~~ 2 decay after a character- 
istic time t c = N 2 t s . The corresponding characteristic 
length scale V Dt c (with diffusion constant D) is the lo- 
calisation length in an iV-mode waveguide. The crossover 
is therefore a dynamical signature of localisation in the 
reflectance of a random medium, distinct from the sig- 
nature in the transmittance (or conductance) considered 
previously in the literature 



Let us first formulate the problem more precisely. We 
consider the reflection of a scalar wave (frequency of) from 
a disordered region (length L, mean free path I = cr s ) 
embedded in an iV-mode waveguide (see Fig. 1, inset). 
We assume that the length L is greater than the locali- 
sation length £ = Nl, so that transmission through the 
disordered region can be ignored. If in addition the ab- 
sorption length is greater than £, the reflection matrix 
r(u>) can be regarded as unitary. The matrix product 

C(u,6u) =r*(uj- \6u) r{uj + \5uj) (1) 

is unitary for unitary r, so that its eigenvalues are phase 
factors exp(i0„). 

The power spectrum for a pulse incident in mode n 
and detected in mode m is related to C by 
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(We have normalised J dta(ui,t) = 1.) The phase shifts 
have the joint distribution function P(<frx, 4>2, ■ ■ . , 4>n)- 
To calculate the average of Tr C it suffices to know the 
one-point function p(4>\) = NJ d<p2 ■ ■ ■ J d<f>N P(4")} since 
(TrC) = / d<P /9(0)e l * . We will present two different 
methods of exact solution. The first method || (based 
on analytic continuation) is simple but restricted to the 
one-point function, while the second method B (based 
on a Fokker-Planck equation) is more complicated but 
gives the entire distribution function. 

Analytic continuation to the imaginary frequency dif- 
ference Sto = i/r a relates exp (i(f) n ) to the reflection eigen- 
value R n of an absorbing medium with absorption time 
r a . The one-point functions are related by 
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R n p(R) dR. 
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This is a quick and easy way to solve the problem, since 
p(R) is known exactly as a series of Laguerre polynomials 



1 



JlQ I . The method of analytic continuation is restricted to 
the one-point function because averages of negative pow- 
ers of exp (i<p n ) are not analytic in the reflection eigen- 
values For example, for the two-point function one 
would need to know the average (exp {i<j> n — i<f> m )) —> 
(RnR^ 1 ) that diverges in the absorbing problem. 

The calculation of the power spectrum from Eqs. (||) 
and (H|) is easiest in the absence of time-reversal sym- 
metry, because p(R) then has a particularly simple form 
|fiO| . We obtain the power spectrum 
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where P n is a Legendre polynomial and 9(t) is the unit 
step function. The coefficient a = 2,7r 2 ,8/3 for dimen- 
sionality d — 1,2,3. In the single- mode case Eq. (^) 
simplifies to 



a(u,t) = 4r t (w)[t + AT s {u)}- 2 9(t), 

which is the result of White et al. . It decays as t 
For N — > oo Eq. (Eh simplifies to 
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o(w,t) = t^expi-t/aTsiw^hit/arsiw)]^), (6) 

where I\ is a modified Bessel function. The power spec- 
trum now decays as t~ 3 / 2 . For any finite N we find a 
crossover from a — y / aT s /2nt~ 3 / 2 for t s « i < N 2 t s to 
a = 2Nar s /T 2 for t > N 2 t s . 

In the presence of time-reversal symmetry the exact 
expression for a(u>, t) is more cumbersome but the asymp- 
totics carries over with minor modifications. In particu- 
lar, the large- N limit (JsJ) with its t -3 / 2 decay remains the 
same, while the t~ 2 decay changes only in the prefactor: 
a = (N + l)ar s t~ 2 for t > N 2 t s . 

We now turn to the second method of solution, based 
on a Fokker-Planck equation for the entire distribution 
function P(c/>). The equation for N — 1 was derived in 
Ref. [§. The multi-mode generalisation can be obtained 
most directly by analytic continuation of the Fokker- 
Planck equation for the probability distribution of the re- 
flection eigenvalues — which is known ||To|] . The resulting 
Fokker-Planck equation for the phase shifts takes a sim- 
ple form in the variable z = In cotan(</>/4) £ (—00,00) 
for 4> S (0, 2n). It reads 
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n(z) = 2_. (In cosh z„ — c/3Asinhz„) 

n=l 

+ (3 In |sinh z m — sinh z n \ . 
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We defined the dimensionless frequency increment A = 
aT s 5u> and abbreviated cp = \(3{N — 1) + 1. The index 
(3 = 1(2) in the presence (absence) of time-reversal sym- 
metry. We emphasise that, although the Fokker-Planck 
equation can be obtained by analytic continuation, its 
solution cannot. Indeed, this would give the solution 
P oc e n , that fails because it is not normalisable. 

We proceed as in Ref. Jll| by substituting P(z) = 
ty(z) exp [if2(,z)] , in order to transform the Fokker- 
Planck equation (Q) into the Schrodinger equation 
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U(z,z>) = b(P-2) ™ Sh2z + C ° Sh2z ' 
2 (sinh z — sinh z') 
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(8b) 



Here V = j^P 2 {N - 1)(N - 2 + 6/(3) + ±. By restricting 
ourselves to (3 = 2, the interaction term U vanishes and 
Eq. (g) has the form of a one-dimensional free-fermion 
problem. The general solution is given by the Slater de- 
terminant 



*(*)= Det{^ n (z m )}", 



(9) 



where ^ s an eigenfunction with eigenvalue [i of the 

single-particle equation 



(10) 



The choice of the eigenfunctions is restricted by the con- 
dition Xm=i ^ n = ® that the total eigenvalue vanishes. 

We are now faced with an impasse: The Schrodinger 
equation ( [To| ) has a real spectrum consisting of bound 
states in the potential well V(z). The bottom of the well 
is positive for sufficiently large A, so that the real spec- 
trum contains only positive /i„'s. How then are we to 
satisfy the condition of zero sum of the eigenvalues? The 
way out of this impasse is to allow for complex eigenval- 
ues. The corresponding eigenfunctions will not be square 
integrable, but that is not a problem as long as the prob- 
ability distribution P{z) remains normalisable. This is a 
new twist to the active field of non-Hcrmitian quantum 
mechanics p"2fl . 

The differential equation (|hJ) is known as the conflu- 
ent Heun equation p3| , but we have found no mention of 
the complex spectrum in the mathematical physics liter- 
ature — perhaps because it was considered unphysical. 
The complex spectrum is constructed by means of a com- 
plete set of polynomials to order N—l, 
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The vector of coefficients g = {gi, g%, . . . , <7at} is an eigen- 
vector with eigenvalue (i of the N x N tri-diagonal matrix 
M. with non-zero elements 
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Since the trace of M is zero, the condition J2 n /^n = is 
automatically satisfied. 

The complex spectrum of Eq. (|l(]) consists of the eigen- 
values (i n with two linearly independent sets of eigenfunc- 
tions, 



ijjjz) — V cosh z exp (— iiVAsinhz)^4 M (sinhz), (13a) 
^"(z) = Vcoshzexp (±NA sinh z)i3 M (sinh z). (13b) 

The functions B^ are related to the polynomials A u by 
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The functions and B^ form a bi-orthogonal set on the 
real axis. We choose the normalisation such that 



dx A^ n (x)B lim (x) = 5„ 



(15) 



The solution if; 1 of the first kind can not be used be- 
cause the resulting distribution P(z) is not normalisable. 
In fact, since 



Det {Af, n (x m )} 
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(16) 



one sees that the substitution of tp 1 into Eq. (JoJ) yields 
the solution P cx e that we had rejected earlier. The 
solution ip 11 of the second kind does give a normalisable 
distribution, 



P(z) cx TT (sinh z n — sinh z m ) J^J cosh z% 

x Det {B^ (sinh z m )}^ m=v (17) 

or, in terms of the variable x = sinhz = cotan (0/2), 

P(*)oc Det{^ n (x m )}^ >m=1 Det{^ fc (xO}£ =1 . (18) 

This is the exact solution of Eq. (0) for j3 = 2. Corre- 
lation functions of arbitrary order can be obtained from 
Eq. (|l8|) in terms of a series of the bi-orthogonal func- 
tions Afj, and B u [|L4 15|. For example, the density of 
eigenphases p(0) is given by 



(19) 



Let us examine this solution more closely in various 
limits. For N — 1 one has (i — 0, A u {x) — constant, and 
we reproduce the known single-mode result [p|,[l6f 

PU) = - Ta Su(l- cos ^Ime^Ei (-(),. (20) 

Here £ = 4iT s &j(l — e 1 ^) -1 and Ei is the exponential- 
integral function. For N > 1 the eigenvalues /i„ re- 
main real for AN 2 <C 1. In this regime the inte- 
gral ( |l4]) is easily evaluated, because one can substi- 
tute effectively [(x — x 1 ) 2 + 1] _1 — > 7T(5(x — x'). Hence 
Bfj,(x) cx exp(-iVAx).4 p (x)6>(x). Using again Eq. (|l|), 
we obtain 
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This is the Laguerre ensemble of random-matrix theory. 
The distribution is dominated by x = cotan (0/2) 3> 1, 
so that one can replace x n — > 2/0 n and recover the result 
fil| that the inverse time delays, l/r„ = lim5 W _>o $w/(f> n , 
are distributed according to the Laguerre ensemble. The 
condition AN 2 <C 1 for Laguerre statistics means that 
the characteristic length Ls u — y/D/du; associated with 
the frequency increment 8lj is greater than the localisa- 
tion length £. We therefore refer to the regime of validity 
of Eq. ( |2"l| ) as the localised regime. 

At the opposite extreme we have the ballistic regime 
L§u < i, or A > 1. The integral ( p"4| ) is now dom- 
inated by x' -C 1, hence B l _ l (x) oc (x + l)~ 1 A~ 1 (x). 
Moreover, the off-diagonal elements of the matrix M 
may be neglected so that the polynomials have a 
simple structure: An n (x) cx (x + i) N ~ n (x — i) n_1 . 
The corresponding functions B are given by B fln (x) = 
sin N+1 (0/2) exp [i<j>(n — N — |)]. The resulting distri- 
bution of the eigenphases in the ballistic regime is 



P(4>) cx J} 



J<P„ 



(22) 



which we recognize as the circular ensemble of random- 
matrix theory |l7j . This is as expected, since for large Slo 
the matrix C is the product of two independent reflection 
matrices r, each of which is uniformly distributed in the 
unitary group. The circular ensemble is the correspond- 
ing distribution of the eigenphases. 

The intermediate regime I <C Ls u <C £, or N~ 2 <C 
A <C 1, is the diffusive one. To study this regime we 
make a WKB approximation of the Schrodinger equa- 
tion (|l0|). This approximation requires N 2 A ^> 1 and 
N 3> 1, hence it contains both the ballistic and the dif- 
fusive regimes. We obtain 

= c( M ) (V(z) - M )- 1/4 e ± ^ d " (23) 
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where c(/i) is a normalisation coefficient. The "+" sign 
in the exponent refers to ip* 1 and the "— " sign to ip*. 
The eigenvalues p n densely fill a curve C in the com- 
plex plane. We may substitute / M — > J c p(p)f(p) dp, 
where p(p) = N [47ric 2 (/z)] _1 is the eigenvalue density. 
For analytic /(//), the integral along C depends only on 
the end points p± = N 2 (^±iA) of the curve. From Eqs. 
(|l3|), ( fl9"| ) and (|23| ) we obtain the eigenphase density 
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This result can be obtained also from Eq. (^) in the limit 
N — > oo. It is derived here for (3 = 2, but is actually (3- 
independent. (The /3-dependent corrections in the diffu- 
sive regime are due to weak localisation, and are smaller 
by a factor 1/N.) One can check that p(0) — ► N/2tt 
for A ^ 1, as expected in the ballistic regime. In the 
opposite regime A 1 it simplifies to 
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TV^/cotan (0/4) 
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with the additional restriction: sin (0/2) > A/4. We 
have plotted Eq. (|24|) in Fig. [I] for several values of A. 




FIG. 1. Density of the eigenphases for different values of 
the dimensionless frequency difference A = ar 3 Suj. The 
solid curves are computed from Eq. (^), the data points 
result from a numerical solution of the wave equation on a 
two-dimensional square lattice (a — ir 2 /4, N — 20; the scat- 
tering time t s was obtained independently from the locali- 
sation length). The inset shows the geometry of a random 
medium (shaded) embedded in a waveguide. 

In conclusion, we have presented a signature of local- 
isation in the decay of the power spectrum of a pulse 
reflected by a disordered waveguide. This result is an 
application of the distribution of the correlator of the re- 
flection matrix at two different frequencies, that we have 
calculated for arbitrary number of modes N, scattering 



time t s , and frequency difference Su> . With increasing 
6u the distribution crosses over from the Laguerre en- 
semble in the localised regime (5u -C 1/N 2 t s ) to the 
circular ensemble in the ballistic regime (5u> 1/t s ), via 
an intermediate "diffusive" regime. The distribution in 
this intermediate regime does not have the form of any 
of the ensembles known from random-matrix theory and 
deserves further study. 

This work grew out of an initial investigation of the 
single-mode case with K. J. H. van Bemmel and P. W. 
Brouwer. We thank H. Schomerus for valuable discus- 
sions. Our research was supported by the Dutch Sci- 
ence Foundation NWO/FOM and by the INTAS grant 
97-1342. 
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